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We present a model of fast hadronization of constituent quark matter in rela¬ 
tivistic heavy ion collisions based on rate equations and capture cross sections 
in non-relativistic potential. We utilize a thermodynamically consistent approach 
with a non-ideal equation of state including correlation terms based on string phe¬ 
nomenology. We investigate strange and non-strange particle ratios observed in 
CERN SPS experiments. 


1 Correlations, Clusters, Confinement 

Quark matter in the moment of hadronization is one of the most strongly 
correlated system we deal with in physics. Confinement forces are long range 
forces, so no ’’ionization” of bound clusters of constituent quarks is possible 
below the deconfinement temperature. In order to work with an ’’effective” 
theory of the strongly interacting quark matter, even in its simplest form of 
constructing an descriptive phenomenological equation of state, we have to 
take into account effects of correlation. 

A comparison with well known hard core problem in nuclear matter may 
be enlightening. Although the underlying pair potential diverges at zero dis¬ 
tance (r —> 0), the relative wave function is zero at this point due to exactly 
this repulsive core. As a consequence the effective pair potential is finite at 
small distance giving rise to finite corrections to the mean field energy (corre¬ 
lation energy). In case of the confining static potential between massive quarks 
the relative wave function vanishes beyond a characteristic distance (the Bohr 
radius of hadronic bound state). Here also the divergence at r —> oo is regu¬ 
larized by the pair correlations and the modification to the mean field energy 
and pressure is finite. 
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The above consideration results in an effective equation of state at the 
characteristic particle density determined by the minimum of the effective pair 
potential. This concept can be, however, extended to all densities: the correla¬ 
tion improved pair potential results in a better high-density behavior of nuclear 
and in a better low-density behavior of quark matter. Eventually a certain, 
nonlinear density-dependence in the equation of state takes such correlations 
into account in a phenomenological way. 

In our present model massive quarks and their hadronic (i.e. color neutral) 
clusters are constituents of a mixture. We describe this mixture with a quasi 
particle Hamiltonian including a density dependent background energy density. 
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These thermodynamical consistency relations are satisfied by the following 
general form of the free energy 

F(iVi, iV2, ...,V,T)=J2 V, T) + E6(iVi/E, N^/V,...). (5) 


Here Ni = Vru is the number of constituents, V the reaction volume and T 
the temperature. We choose the free energy as the thermodynamical poten¬ 
tial, because for a system of ’’quark to hadron transition” - chemical, shortly 
’’transchemical”, equations the time evolution of the numbers can be obtained 
relatively easy, contrary to that of the chemical potentials jii. Note that the 
6(ni, n 2 ,...) correction does not depend on the temperature or on the volume 
explicitly. The derived thermodynamical quantities are 

id dh 


2 



S = S' 


id 


E = E 


lid 


Vb. 


id 7 


( 6 ) 


The relation to the vector mean field Ui and to the background (bag) energy 
density is given by 
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Let us consider a few examples. For the good old MIT Bag model b = B is 
constant, leading toU = B and vanishing vector mean field for the constituents 
Ui = 0. The pressure p = UgT — B = aT'^ — B is negative below a critical 
T = Tc, where the quark matter is unstable. 

In a pure mean field approach b is linear, b = Birii, leading to C/ = 0 
(no bag constant) and Ui = Bi. This system is stable at all temperatures 
(P =P“^)- 

Finally we consider a simple system of massive quarks interacting via a 
string-like potential. In this case the correction to the free energy is propor¬ 
tional to the average distance of quarks, 

b= Kn{£) = . (8) 

In this case one obtains corrections to both mean fields, U = (l/3)(Tn^/^ and 
Uq = (2/3)crn^/^. The correlation correction is distributed between the back¬ 
ground (’’bag” constant) and the vector mean field for the quarks. The pressure 
of this system also becomes negative below a certain temperature: 

p = nT- « aT^ - bT^- (9) 

This property is common with the bag model. The string EOS, however, 
results zero pressure at zero density (so no vacuum renormalization is needed). 
Furthermore - due to the vector mean field correction - the string EOS shows 
a non-ideal behavior of the chemistry, which is absent in the bag model. 

Fig.|^ shows the dependence of the quark chemical potential on the density, 
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Since the entropy production in an isolated expanding system is given by 

S = -J2^N,>0, ( 11 ) 

i 

any physical description of the coalescence rate of quarks possesses the property 
that n and N have opposite signs. In the particular two-body fusion process 
studied in the transchemistry model of next Section, one has 

N oc -iV^ (l - , (12) 

satisfying the above condition. 

Due to this often cited entropy growth principle, the quark density in a 
constant volume decreases whenever /r > 0 and increases whenever /r < 0. 
The isotherms of the ^{n) relation depicted in Fig.|^ reveal that below a given 
temperature Tchem no chemical equilibrium is possible in a pure quark matter 
with string EOS. In this case the quark chemical potential is overall positive at 
any density — so quarks will be eliminated from the system. The only natural 
way to find a thermodynamically stable state is to form color neutral clusters 
which makes no contribution to the string correction in the equation of state. 
This is what happens by the hadronization of colored constituent quark matter 
and this is the basis of the transchemistry model discussed in the next section. 



Figure 1: Isotherms of the chemical potential are plotted against the quark density with 
a non-ideal equation of state with string-like interaction. The infinite temperature case 
coincides with the ideal gas form, while below a critical temperature there is no chemical 
equilibrium possible: the quark density tends towards zero. 


Finally we would like to show some results of a molecular dynamical calcu¬ 
lation of massive quark matter in the framework of the chromodielectric model. 
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Figure 2: Simulation snapshots after initialization of a central-rapidity slice of a QGP tube 
with Lorentz-invariant flow (“Bjorken scenario”) in the chromodielectric model. Colored 
objects (quarks, gluons) are drawn in dark grey, hadrons in light grey. 


This model describes the confining strings due to an abelian Gauss law but 
with a nontrivial dielectric constant. Inspecting the snapshots of quark matter 
evolution in a Bjorken scenario (Fig.^ the ramification is clearly seen in the 
intermediate stage of hadronization. 


2 The transchemistry model 

A detailed description of the transchemistry model can be found in0. The 
equation of state used for the description of the evolution of constituent quark 
matter into hadrons during a sudden expansion utilizes a generalization of the 
above discussed string-like correlation energy. In the hadronizing mixture we 
deal with quarks, anti-quarks, diquarks, anti-diquarks carrying color charge 
and with mesons, baryons and anti-baryons being color neutral. The number 
of strings is assumed to be proportional to a weighted sum of color particles 
Q = We take qi = Q for hadrons, qi = 1 for quarks and qi = 3/2 for 

_ 1 /Q 

diquarks. The average length, {i) = nc is determined by the color density 
ric = excluding hadrons from the sum. The correction to the free 

energy of a mixture of massive ideal gas constituents without inner degrees of 
freedom is then given by 

AF = an-^I^Q. (13) 
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The modification of the non-hadronic chemical potentials becomes 

( 14 ) 

The non-ideal cooling law, equivalent to = 0 for a perfect fluid of the 

mixture, becomes 
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re-heating the system while eliminating color charges. The cooling is due to 
the expansion {V > 0) and due to rest mass production by making hadronic 
resonances. 

The initial temperature and volume was obtained in a one dimensional 
Bjorken scenario from the bombarding energy and the stopping power. At 
the CERN SPS Pb + Pb experiment we obtain Vq « 400 fm^ and Tq w 180 
MeV. The initial quark numbers include direct and gluon fragmentation pair 
creation and sum up to about 3000 with a strangeness ratio of fs = 0.21. 

The general two body fusion process is supported by the in-medium cross 
section of massive constituents in a Coulomb potential in medium. This is 
proportional with a screening volume, p^, which also depends on the presence 
of color charges. We make the assumption that p^ric is constant obtaining a 
dynamical confinement effect of increasing hadronization cross sections in the 
low color density stage. These phenomenologically amounts to the string-like 
correlation correction to the hadronization cross sections and influences the 
dynamics of the transchemical evolution. 

Each chemical reaction is of type i -I- j fc, accounting for the changes 
dNi = dNj = —Adt and dNk = Adt in the particle numbers. The general 
form, 

A = (l-e"^-^-^) , (17) 

assures the increase of entropy while approaching chemical equilibrium. Fi¬ 
nally, in order to get the final hadron numbers, hadronic decays are taken 
into ^count with the dominant branching ratios obtained from Particle Data 
Tablfl. 
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3 Constituent Quark Matter at SPS 

In this section we present results of the numerical solution of 43 coupled differ¬ 
ential equations describing the chemical and thermal evolution of the quark, 
diquark - hadron mixture. In FigJ| the evolution of the temperature, entropy, 
pressure and energy density is shown for a 168 GeV/nucleon Pb + Pb collision. 
At the beginning there is a rapid drop of the temperature (cf. Fig.||a) due 
to heavy hadron resonance formation. It is followed by a mild re-heating as 
an effect of the color confinement (by producing a colorless hadron the kinetic 
energy of the relative motion and the string energy is removed and replaced by 
the rest mass of the hadron). The total entropy monotonically increases during 
the hadronization while the contribution of colored particles is gradually elim¬ 
inated (FigJ|b). The partial pressure of the constituent quark plasma rapidly 
decreases, it becomes even negative as the color density drops (Fig.^). The 
total pressure, however, remains positive - so the mixture is stable and contin¬ 
ues the adiabatic expansion. The partial and total energy densities evolve in 
the expected way going over to the e oc 1/r behavior at very late times only, 
when p = 0. 




Figure 3: Evolution of thermodynamical quantities in the transchemistry model 

Fig.^ presents the time evolution of constituent particle numbers: quarks, 
diquarks (left column from top to bottom), mesons and baryons (right column). 
The evolution of anti-particles is qualitatively similar, their absolute numbers 
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are only smaller. The hadron numbers shown here are those before hadronic 
resonance decays. It is clear from these figures that mesons are formed faster 
than baryons, the reason being that in this model the baryons are produced by 
a two step process and the initial diquark content was zero in these particular 
calculations. 


Table 1: hadron multiplicities for Pb+Pb collision at 168 GeV/nucleon bombarding energy. 


Pb+Pb 

NA49 

TrCHEM. 

ALCOR 

RQMD 

K+ 

76 

76.14 

78.06 

79.0 

K- 

32 

38.36 

34.66 

50.4 


54 

57.25 

56.36 

63.5 

P-P 

145 

151.6 

147.03 

171.8 

AO-like 

50+ 10 

62.30 

69.07 

56.8 

AO-like 

8 ± 1.5 

8.14 

8.12 

19.3 


In the ALCOR model 0 the ratio of hadronic species are determined by 
the steepness of these curves (i.e. by Ni{0)). Since the Ni{t) curves do not 
cross, the algebraic approach of ALCOR to the solution of the underlying rate 
equations works well. 




Figure 4: Evolution of thermodynamical quantities in the transchemistry model 


Finally in Table ^ some hadron numbers obtained in the transchemistry 
model are compared with experimental data of the NA49 collaboration and 




































with two other theoretical models: 
comparison can be found inld. 


ALGOR and 


RQMri. A 


more detailed 


4 Conclusion 

We presented a new model for the hadronization of constituent quark plasma 
based on a string-like correlation energy in the equation of state and on rate 
equations in a quark matter - hadron matter mixture. Due to large hadroniza¬ 
tion rates and due to a very off-equilibrium initial state, over-saturated with 
(massive) quarks, the hadronization process is fast: its time-scale is a few 
fm/c. We also found that in the particular situation of 168 GeV/nucleon 
heavy-ion collisions the algebraic ALGOR model is a good approximation to 
the numerical solution of the transchemical rate equations. The comparison 
with presently available experimental data indicate, that in such collisions a 
piece of matter is formed, inside which massive quarks, diquarks and their anti 
particles interact with a string-like mean field. 
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